Detecting gravitational waves from precessing binaries of spinning compact objects. 
II. Search implementation for low-mass binaries 
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Detection template families (DTFs) are built to capture the essential features of true gravita- 
tional waveforms using a small set of phenomenological waveform parameters. Buonanno, Chen, 
and Vallisneri [Phys. Rev. D 67, 104025 (2003)] proposed the "BCV2" DTP to perform computa- 
tionally efficient searches for signals from precessing binaries of compact stellar objects. Here we 
test the signal-matching performance of the BCV2 DTP for asymmetric-mass-ratio binaries, and 
specifically for double-black-hole binaries with component masses (mi,m2) € [6, 12]M0 x [1,3]M0, 
and for black-hole-neutron-star binaries with component masses (mi,m2) = (IOMq, 1.4Mq); we 
take all black holes to be maximally spinning. We find a satisfactory signal-matching performance, 
with fitting factors averaging between 0.94 and 0.98. We also scope out the region of BCV2 pa- 
rameters needed for a template-based search, we evaluate the template match metric, we discuss a 
template-placement strategy, and we estimate the number of templates needed for searches at the 
LIGO design sensitivity. In addition, after gaining more insight in the dynamics of spin-orbit pre- 
cession, we propose a modification of the BCV2 DTP that is parametrized by physical (rather than 
phenomenological) parameters. We test this modified "BCV2P" DTP for the {IQMq,1AMq) black- 
hole-neutron-star system, finding a signal-matching performance comparable to the BCV2 DTP, 
and a reliable parameter-estimation capability for target-binary quantities such as the chirp mass 
and the opening angle (the angle between the black-hole spin and the orbital angular momentum). 

PACS numbers: 04.30.Db, 04.25. Nx, 04.80.Nn, 95.55.Ym 



I. INTRODUCTION 

As ground-based gravitational-wave (GW) detectors 
based on laser interferometry approach their design 
sensitivities, the emphasis in data analysis is shifting 
from upper-limit studies to proper detection searches. 
In addition, the length of data-taking runs is stretching 
to several months, with typical duty cycles approaching 
unity, substantiating the need for online (or at least real- 
time) searches to be performed as data become available. 
It is then crucial to develop search algorithms that maxi- 
mize the number of detections while making efficient use 
of computational resources. 

Inspiraling binaries of black holes (BHs) and/or neu- 
tron stars (NSs) are among the most promising and 
best-understood sources for GW interferometers, which 
can observe the waveforms emitted during the adia- 
batic phase of these inspirals, well described by post- 
Newtonian (PN) calculations [J. For these signals, the 
search algorithms of choice are based on matched filter- 
ing 0, whereby the detector output is compared (i.e., 
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correlated, after noise weighting) with a hank of theoret- 
ically derived signal templates, which encompass the GW 
signals expected from systems with a prescribed range of 
physical parameters. 

Reference introduced the phrase "detection tem- 
plate families" (DTFs) to denote families of signals that 
capture the essential features of the true waveforms, but 
depend on a smaller number of parameters, either phys- 
ical or phenomenological (i.e., describing the waveforms 
rather than the sources). At their best, DTFs can re- 
duce computational requirements while achieving essen- 
tially the same detection performance as exact templates; 
however, they are less adequate for upper-limit studies, 
because they may include nonphysical signal shapes that 
result in increased noise-induced triggers, and for param- 
eter estimation, because the mapping between template 
and binary parameters may not be one-to-one, or may 
magnify errors. In Ref. 0, the "BCVl" DTP was de- 
signed to span the families of nominally exact (but par- 
tially inconsistent inspiraling-binary waveforms obtained 
using different resummation schemes to integrate the PN 
equations. 

A reduction in the number of waveform parameters is 
especially necessary when the binary components carry 
significant spins not aligned with the orbital angular mo- 
mentum; spin-orbit and spin-spin couplings can then in- 
duce a strong precession of the orbital plane, and there- 
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fore a substantial modulation of GW amplitude and 
phase 0- Detection-efficient search templates must ac- 
count for these effects of spin, but a straightforward 
parametrization of search templates by the physical pa- 
rameters that can affect precession results in intractably 
large template banks. 

To solve this problem, several DTFs for precessing 
compact binaries have been proposed in the past decade 
[llElSIISIilliilll- A DTF based on the Apostolatos 
ansatz for the evolution of precession frequency was am- 
ply investigated in Refs. 11. 13.j , and an improved version 
was proposed in Ref. However, the computational 

resources required by the Apostolatos-type families are 
still prohibitive; more important, their signal matching 
performance (i.e., the fitting-factor FF) is not very sat- 
isfactory. 

In Ref. 01, Buonanno, Chen, and Vallisneri analyzed 
the physics of spinning-binary precession and waveform 
generation, and showed that the modulational effects can 
be isolated in the evolution of the GW polarization ten- 
sors, which are combined with the detector's antenna pat- 
terns to yield its response. As a result, the response of 
the detector can be written as the product of a carrier 
signal and a complex modulation factor; the latter can be 
viewed an extension of the Apostolatos formula. In Ref. 
[T^ . the precessing waveforms were cast into a math- 
ematical form (the linear combination of three simpler 
signals, with complex coefficients) that allows searching 
automatically and economically over all the precession- 
related parameters, except for a single parameter B that 
describes the timescale of modulation. Henceforth, we 
shall refer to the template family proposed in Ref. |lj| as 
the "BCV2" DTF. 

In Ref. [13, the BCV2 DTF was tested for precessing 
BH-BH binaries with high total mass (12Afo < M < 
30AfQ) and comparable component masses, and for the 
single mass configuration (10+ \A)Mq, representative of 
NS-BH systems. In all cases, the signal-matching per- 
formance was good (FF > 0.9), with consistent improve- 
ments over search templates that do not include preccs- 
sional effects (for instance, in the NS-BH system the FF 
increases from ~ 0.78 to ~ 0.93). Signals from precessing 
binaries with asymmetric component masses are harder 
to match, because they have more orbital and preces- 
sional cycles (i.e., more complex waveforms) in the band 
of good interferometer sensitivity. 

In this paper, we extend the BCV2 performance analy- 
sis of Ref. [U to asymmetric mass ratios, taking into con- 
sideration systems with component masses (7711,7712) G 
[6, 12]M0 X [1,3]Mq, for which we expect a large num- 
ber of precession cycles (see Fig. below). In addition, 
we estimate the region of the DTF parameter space that 
must be included in a search for such systems; we calcu- 
late the template match metric [ij, [T^ [13I ; we provide 
a strategy for template placement; last, we estimate the 
number of templates required for the search. After recon- 
sidering the Apostolatos ansatz, we are also able to shed 
new light on the phcnomcnological parameter B that de- 



scribes the timescale of modulation; indeed, we derive an 
explicit formula for the evolution of the precession angle 
in terms of the physical parameters of the binary, and we 
use this formula to propose a modification of the BCV2 
DTF that dispenses with B. 

While this paper is concerned with DTFs for precess- 
ing binaries, we note that a physical template family for 
single-spin precessing compact binaries wa,s proposed in 
Ref. [Tjl, and thoroughly tested in Ref. [ig. The at- 
tribute "physical" is warranted because the family is ob- 
tained by integrating the FN equations Q in the time 
domain, and the templates are labeled by the physi- 
cal parameters of the binary. Furthermore, Ref. [T3 
showed that the single-spin physical family has a satisfac- 
tory signal-matching performance also for the waveforms 
emitted by double-spin precessing compact binaries, at 
least for component masses (7711,7712) € [3,15]AfQ x 
[3, ISJM©; moreover, the parameters of the best-fit single- 
spin templates can be used to estimate the parameters of 
the double-spin target systems [Tgl l. However, this phys- 
ical template family may be more complicated to imple- 
ment and more computationally expensive (and there- 
fore less attractive for use in online searches) than the 
frequency-domain DTFs such as BCV2. 

This paper is organized as follows. In Sec.m we briefly 
review the BCV2 DTF and the Apostolatos ansatz, and 
we discuss how the phenomenological parameter S, which 
describes the timescale of precession, can be related to 
the physical parameters of the binary. In Sec. IIII Al we 
discuss the signal-matching performance of the BCV2 
DTF over a range of binary component masses. In 
Sec. IIII B1 we introduce a version of the BCV2 DTF mod- 
ified to include the physical evolution of the precession 
angle in single-spin binaries, and we test its performance 
for NS-BH inspirals. In Sec. IIV Al we compute the tem- 
plate match metric for the BCV2 DTF. In Sec. HVBl we 
provide a strategy for template placement, and we esti- 
mate the number of templates required in a search. Last, 
in Sec. [V] we summarize our conclusions. 

In the following; the binary component masses are de- 
noted by 7?7i and 777,2 (with 7?7i > 7772); the symmet- 
ric mass ratio and the total mass by 77 = mim2/M'^ 
and AI = mi + 7)72; the binary component spins by 
Si = Xi ™i E^iid S2 = X2 "73. For single-spin binaries, we 
assume Si = X7?7i and S2 = 0. Throughout the paper, 
the signal-matching performance of DTFs is evaluated 
against a target model for processing binaries governed 
by Eqs. (6)-(32) of Ref . ^3; this target model is valid in 
the adiabatic phase of inspiral, when dynamics arc cor- 
rectly described by FN equations. We use an analytic 
fit to the LIGO-I design noise spectrum (given, e.g., by 
Eq. (28) of Ref. 0); we adopt the standard formahsm 
of matched- filtering GW detection; we follow the con- 
ventions of Ref. {ljj| , which contains a useful glossary of 
matched-filtering notions and quantities; last, we always 
set G = c = 1. 
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II. FEATURES OF PRECESSION DYNAMICS 
IN SINGLE-SPIN BINARIES 

A. Review of the Apostolatos ansatz and of the 
BCV2 DTF 

Apostolatos, Cutler, Sussman, and Thorne (ACST) 
investigated orbital precession in binaries of spinning 
compact objects in two special cases: (i) equal-mass bina- 
ries (toi = 7712), where the spin-spin coupling is switched 
off, and (ii) single-spin binaries (5*2 = 0). In these cases, 
precessional dynamics can always be categorized as sim- 
ple precession or transitional precession. In simple pre- 
cession, the direction of the total angular momentum J 
is roughly constant, while the orbital angular momen- 
tum L and the total spin S = Si -I- S2 precess around it. 
ACST were able to derive an analytical solution for the 
evolution of simple precession (see Sec. IV of Ref. 0]). 
Transitional precession occurs when, during evolution, 
L and S have roughly the same amplitude and become 
nearly antialigned. When this happens, |J| is almost zero 
and J can change suddenly and dramatically. Although 
transitional precession is too complicated for analytical 
treatment, it occurs rarely 0.ll4j. so we will ignore it in 
this paper. 

GW signals from generic precessing binaries are well 
approximated by simple-precession waveforms when the 
ACST assumptions are valid, which happens for two 
classes of binaries: (i) BH-BH binaries with compara- 
ble component masses where the spin-spin interaction 
can be neglected, which are equivalent to systems where 
a single object carries the total spin of the system; (ii) 
BH-NS or BH-BH binaries with very asymmetric mass 
ratios, which can be approximated as single-spin systems 
because the spin of the lighter object is necessarily small. 
It is not guaranteed a priori that simple-precession wave- 
forms can describe also signals emitted by BH-BH bi- 
naries with intermediate mass ratios and/or important 
spin-spin effects. However, it was recently shown [T9| 
that simple-precession waveforms are adequate also for 
these classes of binaries, although the dynamical evolu- 
tion of L and S can exhibit rather different features [3 . 

In simple precession, L and S precess around J (which 
is roughly fixed) with the precession frequency given by 
Eq. (42) of Ref. 0, 



40 



30 

C 



K=0.9 



20 



10 







50 



K=0 



K=-0.9 

100 150 200 250 300 350 400 450 
Frequency (Hz) 



FIG. 1: Evolution of precession angle (radians) as a function 
of GW frequency, for (10 + 1.4)M0 binaries with opening 
angles k — 0.9, and —0.9. The dotted curves show the nu- 
merical evolution of our target systems, while the continuous 
curves follow the analytical expression 



where lu is the orbital angular frequency, and r is the 
orbital separation, so it is straightforward to get ap = 
Bif^^; this regime corresponds to comparable-mass bi- 
naries, or to binaries at large separations (i.e., in the 
early stages of the inspiral). When L S and J ~ 5, 
we have ap = B2f~^^^; this regime corresponds to bina- 
ries with large mass asymmetry, or to binaries at small 
separations (i.e., in the late stages of the inspiral). The 
analytical expressions for the coefficients Bi and B2 are 
given by Eq. (45) of Ref. 0, and depend only on the 
masses and on the total spin of the binary, but not on the 
opening angle between the spin and the orbital angular 
momentum. Although the power laws were derived for 
simple precession under the ACST assumptions, it turns 
out that they can model the dynamics of more general 
configurations, as shown in Refs. [T^H^ . 

On the basis of the ACST analysis, Apostolatos 
reasoned that GWs from precessing binaries should be 
modulated by the orbital precession frequency f2p, and 
suggested adding an fip-depcndent modulation term to 
the nonspinning waveform phasing. 
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ACST found that the evolution of the precession angle 
Up can be approximated by power laws in two extreme 
cases. When L ^ S and J ci: L, using leading-order 
Newtonian expressions, we find 



L^r,M\Mu:)-^'\ r= (^^^ ' , 
0}-^ 5 



(1) (this is the Apostolatos ansatz). Although the result- 
ing DTF has higher fitting factors with precessing-binary 
waveforms than nonspinning DTFs, it is not completely 
satisfactory , especially because of the huge com- 
putational cost implicit in adding the three parameters 
C, 5, B, which are intrinsic {i.e., they increase the dimen- 
sionality of search template banks) . 

In Ref. ^3 J Buonanno, Chen, and Vallisneri proposed 

(2) a DTF (BCV2) based on a modification of Eq. @, 
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^(V'0,'/'3/2,S,/cut,Cfc;/) = 



/ 



-7/6 



(Ci + 1C2) + (C3 + iCi) cos(S/-2/3) + (C5 + ^C6) siii(S/-2/3)l - /)e2-/*o X e'[^o/-=/^+^3/2/-=^^] . (4) 



here to is the signal's time of arrival, and /cut is the cutoff 
frequency; precessional effects are modeled by the modu- 
lation terms cos(S/~^/'^) and sin(S/^^/'^) in the complex 
amplitude, separately from the nonspinning evolution of 
phase. Possible modulation morphologies are enriched 
by the presence of the complex linear-combination coef- 
ficients C3 -|- zC4 and C5 + iC^, improving the efficiency of 
matching to target waveforms. Indeed the complex mod- 
ulation terms modulate both amplitude and phase of the 
nonspinning signal. An important feature of the BCV2 
DTF is that little computational cost is added by the 
search over the parameters Ci^...^6- which are extrinsic 
{i.e., the detection statistic can be maximized analyti- 
cally over them). 

Although the /"^/'^ power law adopted in the BCV2 
DTF Q is expected to be vaUd only for S* > L, (i.e., 
for binaries with asymmetric mass ratios, or in the late 
stages of inspiral) , the DTF yields high fitting factors also 
for comparable-mass binaries, as verified in Ref. [1^ by 
Monte Carlo simulation. The reason is probably to be 
found in the broad variety of modulation morphologies 
parametrized by the B and Ci_..._6 parameters. 



except in the limits L ^ S and L S, where the power 
laws Up = Bif~^'~^^^ do not include k. Let us see what 
function we should expect. From Eq. (49) of 0, we have 
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which was obtained by expressing the total angular mo- 
mentum J in Eq. in terms of the orbital angular mo- 
mentum L. In Eq. |SJ), j{t) denotes the quantity S/L{t); 
the dependence of flp on k vanishes with 7 3> 1 (i.e., 
with L < S") or 7 < 1 (i.e., with L > S"). Using the 
leading-order Newtonian expressions for L and r given in 
Eq. we can integrate flp analytically and obtain ap 
as a function of mi, m2, k, and X- 
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5 4toi + 3m2 



384 TOi 
[~A [(2-3k2)xm + kXa/w"'+2i;-2] 

+ 3k (1 - K^) xii log [kxm + + a] I + const. (6) 

where 



B. Analysis of the DTF parameter B 

In this section we relate the phenomenological preces- 
sion parameter B to the physical parameters of the binary 
(the two masses, 7711,7712, the cosine between the direc- 
tions of the total spin and the orbital angular momentum, 
K = L • S, and the magnitude of the spin, x = Si/ml). 
In doing so, we clarify why the BCV2 DTF is capable of 
mimicking the precessional effects in the target signal. 

In Fig. 16 of Ref. 0|, the distribution of the best-fit 
DTF B is plotted against the target-system parameter 
K. The target waveforms were generated at 2PN order, 
for BH-NS binaries of component masses (10 -I- 1.4)Mq 
with maximally spinning BHs. The spread of the data 
points corresponds to uniform distributions of the initial 
spin and angular-momentum orientations. The points 
seem to cluster around three lines, but no explanation 
is offered for this interesting feature. We are now able 
to explain this behavior; what we learn in the process 
will enable us to construct an improved DTF (BCV2P) 
parametrized by physical parameters (see Sec. IIII 51 . 

Although the best-fit DTF parameter B is not, strictly 
speaking, physical, it is clearly related to the evolution of 
the precession angle Up in the target system. Moreover, 
we expect the best-fit ;B to be a function of the target- 
system opening angle k (as seen in Fig. 16 of Ref. 
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and where / = uj/tt. The 'W in stands for "Newto- 
nian" . These expressions are equivalent to Eqs. (63a) and 
(63b) of Ref. 0- Note that an analytical expression in- 
cluding higher PN corrections could also be given. How- 
ever, for simplicity we prefer to restrict ourselves to the 
lowest order. In Fig.^ for a binary of mass (1O-|-1.4)M0 
and for several values of k, we compare the analytical 
precession angles (/) with numerical values obtained 

from our target models by projecting L(/) onto the plane 
perpendicular to the vector J (which is constant because 
we only consider simple precession), and recording the 
cumulative angle swept by the projected image. We see 
that our leading-order formula reproduces the shape of 
the numerically-obtained curve, although the quantita- 
tive difference is appreciable. This is due to the fact 
that we write at the Newtonian order; the agreement 
would otherwise be perfect. 

We can now try to explain the dependence of the DTF 
parameter B on the parameters of the binary, and the 
clustering seen in Fig. 16 of Ref. j^^. Since the power 
laws Bf~^^^'^ cannot match ap{f) exactly, we estab- 
lish a correspondence by requiring that the instantaneous 
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FIG. 2: Best-fit values of the BCV2 DTF parameter B 
as a function of tlie target-system opening angle k, for 
(1O-|-1.4)M0 binaries with uniformly distributed initial orien- 
tations of spin and angular momentum. The larger connected 
dots show the equivalent B as evaluated from Eq. @ using 
the numerical ap{f) and the analytical Op (/), respectively. 



rate of change of the two precession angles be equal at 
the approximate frequency of best detector sensitivity, 
~ 150 Hz (appropriate for initial LIGO): 

|sr2/3 ^ f7* = ^ , /= 150 Hz, (8) 

and therefore 



and Fig. 16 of Ref. ^^-^ The linear cluster identified 
vifith rip = ilp in Fig. |21 corresponds in fact to the top 
linear cluster in Fig. 16 of Ref. and the fraction 

of points in this cluster is significantly larger here than 
there. Moreover, the top linear cluster corresponding to 
Hp = 2 ftp cannot be seen clearly in Fig. 16 of Ref. [lj|. 
These differences are due to the adoption, for this pa- 
per, of a better numerical overlap-maximization scheme, 
which has a better chance of finding the true global maxi- 
mum overlap. We shall discuss the maximization scheme 
in detail in Sec. IhTH ] 



C. Higher harmonics in templates and signals 

We shall now consider why multiple clusters appear in 
Fig. 13 Wc shall see that the precession frequencies Hp 
and Hp defined in Eq. (jHJ are not usually the only mod- 
ulation frequencies to appear in the spectra of the tem- 
plate and signal waveforms — indeed, they may not even 
be the dominant frequencies. It is then conceivable that 
if different harmonics of the precession frequency domi- 
nate in the template and signal waveforms, the maximum 
overlaps may occur when 17* and Qp arc not equal, but 
instead related by integer factors. 

We first consider the frequency content of the BCV2 
template phase modulations. As discussed above, their 
precession angle is ap = Bf"^^^, and the precession fre- 
quency flp is given by the time derivative of ap. Modu- 
lation effects arc included by way of the complex factor 



6: 



3 .5/3 dapif) 

r df 



f=150Hz 



Here ftp denotes the template's equivalent precession fre- 
quency, while Hp denotes the target's orbital precession 
frequency at 150 Hz. Figure |21 shows a test of the cor- 
respondence, performed for (10 + 1.4)Mq binaries with 
uniformly distributed values of n and of the initial ori- 
entations of total spin and orbital angular momentum. 
The small dots indicate the best-fit values of B, simi- 
larly to Fig. 16 of Ref. The larger dots mark the 
pairs (k, B{k)) obtained from Eq. using the analytical 
ap (/) and the numerical ap{f), respectively. The pairs 
line up quite well with the linear cluster in the middle of 
Fig. 121 which includes the majority of data points. Thus, 
it is correct to state that the value of B represents the 
rate of change fi* = Hp of the precession angle in the 
middle of the frequency band of good detectory sensitiv- 
ity. Similar conclusions hold also for binaries with masses 
{mi, 1712) G [6, 12]Mq X [l,3]Af0. The other linear clus- 
ters seen in Fig. |21 correspond roughly to 2Hp ~ Hp and 
Hp ~ 2Hp, for reasons explained in the next section (the 
clusters in the lower corners, on the other hand, corre- 
spond to systems where the effects of precession on the 
waveforms are negligible). 

[Note that there arc some differences between Fig. |2 



(9) (Ci + iCa) + {C3 + iCi) cos(6/-2/3) + {C^ + iCq) sin(6/-2/3) 

(10) 

(with Ci_..._6 G lR)j which is clearly a periodic function 
of Up , but is obviously far from a simple sinusoid with a 
single frequency. Recasting the factor as 



iCi+iC2)+Ar cos(6/ 



-2/3 



a, C0s(6/-2/3 -f <^,) 
(11) 

(with Ci^2,^r,i, fr.i G K), WC scc that it traces an ellipse 
in the complex plane as ap varies from to 2tt. The 
shape of the ellipse is determined by Ar,i and fr,i, and 
the displacement of its center from the origin by Ci +iC2- 
For two choices of the modulation parameters Ci^2, 
Ar,i, and ipr;i, Fig. 121 shows the complex-plane trajec- 
tory of the modulation factor (left panels), the oscilla- 
tory part^ of the phase as a function of ap = Bf~^^^ 
(middle panels), and the amplitude of its Fourier trans- 
form. In the first example, we choose Ci^2 to dominate, 



^ In Fig. 16 of |14| . the x axis refers to Keff/'^JJff > where K^g = 
L-Scff/M^, with Seff defined by Eq. (7) of Ref. In the limit 
in which only one of the two bodies carries spin, the quantity k 
used in this paper equals K^ff/n^'^. 

^ That is, the residual obtained after fitting the total phase to the 
nonspinning phase 4>o + ^Trtof + '4>of~^^^ + '03/2/~^^^- 
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FIG. 3: BCV2 DTF complex modulation factor, as given by Eq. lllll . The top (bottom) row refers to a Ci,2-dominated 
(^r,i-dominated) choice of the modulation parameters. The left column shows the complex-plane trajectory of the modulation 
factor when the precession angle ap varies between and 27r; the middle column shows the oscillatory part of phase as function 
of Qp; the right column shows the log-amplitude of the Fourier transform of phase. 



setting Ci 1, C2 = 0.5, Ar = 0.2, A., = 0.1, ipr = 0, and 
'fi = 7r/2. The phase is periodic in cip with period 27r, 
and it is rather close to a single sinusoid (upper middle 
panel). The Fourier spectrum (upper right panel) shows 
that the dominant template modulation frequency is 
the precession frequency f2* , and that the contributions 
from higher harmonics are at least an order of magnitude 
smaller. 

In the second example, we choose Ar,i to dominate, 
setting Ci = 0, C2 = 0, Ar = 1, Ai ~ 1.5, Lpr ~ 0, and 
ipi = 7r/2. In this case the modulation ellipse encloses 
the origin, so the phase contains a monotonic component 
that grows by 27r each time the modulation factor com- 
pletes an orbit around the origin. It is easy to see that 
this monotonic component is simply ap. The oscillatory 
part of the phase can be obtained from Fig. O (lower 
left panel) by taking the phase difference between points 
with the same real parts on the ellipse and on the unit 
circle. In this case, the phase is (almost) periodic in ap 
with period tt, as shown in the lower middle panel. The 
Fourier spectrum (lower right panel) shows clearly that 
the dominant frequency is fi^ = 2ri* , and that there is 
no component at 17^ or at any other odd harmonic. 

Investigating other choices of the modulation parame- 
ters, we find that the Fourier spectra of the phase mod- 
ulations have always their highest peaks at either fi* or 
2f2*, depending on whether the modulation ellipse en- 
closes the origin or not. Most of the differences between 
the various cases lie in the structure of harmonics above 
the second. For instance, if the trajectory of the complex 
amplitude has very large ellipticity, higher-order harmon- 
ics can become comparable to the lowest harmonic (i.e.. 



the fundamental or the second harmonic). 

We now consider the frequency content of the target- 
waveform phase modulations. Under the precession con- 
vention introduced in Ref . |l4j , the frequency-domain ex- 
pression for the waveforms, at the leading order in the 
stationary-phase approximation, are given by Eq. (83) of 
Ref. or namely 

h..s^{f) = -hcif) {Mtf)t+^[<^Atf)f} 

X ([T+(i;)]^.,F+ + [Tx(t/)]^.,i^x) for/>0, (12) 

where hc{f) is the unmodulated carrier signal, the 
T-|_.x(i/) are the detector polarization tensors in the ra- 
diation frame [see Eq.(25) of BCV2], and all precessional 
effects are isolated in the evolving GW polarization ten- 
sors e_|__x(*/)- These are defined as 

e+ = ei (g) ei - 62 ® e2 , 

ex = ei (g) 62 + 62 (8) ei , (13) 

where ei and 62 form a basis in the instantaneous orbital 
plane. The time dependence of these polarization tensors 
enters the waveform through terms of the form 

[e+,x(t/)]'MT+,x(t/)],,. , (14) 

which can be approximated by C+,x cos{Bf'^^^ + (5+,x), 
adopting the Apostolatos ansatz. However, this is only 
an approximation, in two distinct ways. First, the com- 
ponents of the polarization tensors depend quadratically 
on the components of the basis vectors ei.2; as a conse- 
quence, even if the basis vectors contained only oscilla- 
tions with frequency fi* the waveform would still end up 
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FIG. 4: Oscillatory part of phase modulation in target waveforms for BH-NS binaries with (mi,m2,x) = {WMq, IAMq, 1), 
and with directional parameters such that the detector line of sight is perpendicular to the initial L-S plane and to the detector 
plane, and that the detector is oriented along the "+" GW polarization. The three panels correspond to k = 0.866 (for which 
L is close to CIl), k, = 0.050, and k = —0.803 (for which L is perpendicular to flz, at / = 30 Hz). The sharp turns observed in 
the second and third panels are not singularities of the phase, but happen when the projections of e+ and ex on the detector 
frame (and therefore the real and imaginary components in Eq. 1121 ') are both small, so that the phase can change very rapidly. 



with both Op and 2f2p modulations. Second, under the 
precession convention, the basis vectors ei.2 [and thus 
the quantities in Eq. (|14|) ] do not really precess together 
with the orbital plane around J with angular velocity 
= Op J; instead, they precess around the component 
of ri/^ that is orthogonal to the orbital angular momen- 
tum L [Eq.(72) of BCV2], 



(15) 



which in turn precesses together with \^tq(t) around 
Q,L{t)- For this reason, the oscillations in ei_2 are 
more complicated than simple sinusoids of frequency Op. 
ACST had already observed this fact, noticing that the 
phase modulations are certainly due to the precessional 
evolution of the orbital plane, but arise also from the so- 
called Thomas precession term (see Eq. (29) of Ref. j3 
and the discussion around it). 

In Fig. 21 we show the modulation of the target- 
waveform phase as a functions of ap for three binary con- 
figurations. As for the discussion of template waveforms, 
we consider only the oscillatory part of phase modula- 
tion. This is done here by fitting the total phase to a 
nonspinning phase + 2ntof + V'o/"^^"^ + ^'3/2/"^^'^: 
and taking the residual. In the first example on the left 
of Fig. 01 the phase modulation is periodic in ap with pe- 
riod 27r, and it is not very different from a single sinusoid. 
So in this case the dominant modulation frequency is the 
precession frequency Op. In the second example in the 
middle of Fig. 01 the dominant frequency is still Op, but 
the 20p component is also quite significant. In the third 
example on the right of Fig. 01 the dominant frequency is 
40p. Other systems with values of k between those used 
for Fig. 01 show similar features, sometimes with a larger 
number of frequency components. 

So far, this discussion suggests that the phase- 
modulation frequencies O*^ and O^^ of the template and 
target waveforms are not always the precession frequen- 
cies Op and Op (which depend on B and on the phys- 



ical parameters of the target system, respectively), but 
can also be their integer multiples. In light of this, how 
can we understand the multiple clustering seen in Fig. ^ 
The answer is that, for a given target signal, there can 
be several templates whose overlap with the target is a 
local maximum, corresponding to different combinations 
of the Op and Op harmonics. 

Suppose for example that the target waveform con- 
tains the first and second harmonics of the precession 
frequency. Op and 20p; then a template with Op = Op/2 
could match these two components with its second and 
fourth harmonics, a template with Op = Op could match 
them with its first and second harmonics, while a tem- 
plate with Op = 20p could match only the Of„ = 20p 
component with its first harmonic. Using this reasoning, 
we can easily understand the existence of clusters of local 
maxima with O^ = O;, 20^, 30^, . . . , and O^ = 0^/2, 
Op/3, Op/4, . . . However, if the DTF is able to reproduce 
the entire harmonic structure of the signal, then the lo- 
cal maximum with Op = Op must also be the global 
maximum. The fact that in Fig. |21 we have clusters at 
O^ = 0^/2 and O* = 20^ suggests that the BCV2 DTF 
cannot do this perfectly. The analysis of the best-fit 
template modulation parameters for the points in the 
Op = Op/2 cluster of Fig. |21 (extending from k ~ —0.2 
to K ~ 1.0) confirms the template frequency-doubling 
scenario discussed in this section. 



III. SIGNAL-MATCHING PERFORMANCE OF 
THE BCV2 AND BCV2P DTFS 

A. Performance of the BCV2 detection template 
family 

The basic diagnostic of DTF signal-matching perfor- 
mance is the fitting factor FF (0 < FF < 1), defined as 
the match between a given template in the target family 
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FIG. 5: Distribution of BCV2 DTF fitting factors for populations of target systems with component masses (7711,7722) = 
{6, 8, 10, 12}Mq X {1, 2, 3}A'f0, and with uniformly sampled directional and local angular parameters. For each pair of masses, 
we include 200 target systems. The curves show the number of samples falling within each bin marked on the abscissa; the bins 
have equal width, but are plotted logarithmically to emphasize FFs close to unity. The figures show also the average fitting 
factor FF and the effective average fitting factor FF^ff (as defined by Eq. (119) of Ref. [Tl|'). with their estimated Monte Carlo 
error in parentheses. 



and the templates in the DTF, maximized over all the pa- 
rameters of the DTF. The maximization of the match in- 
duces a many-to-one projection between the target-signal 
parameters and the best-fitting DTF parameters. The 
way in which the maximization is carried out informs 
the distinction between intrinsic and extrinsic DTF pa- 
rameters: full templates must be recomputed for each 
value of the intrinsic parameters to be explored, while 
the maximum over the extrinsic parameters can be com- 
puted analytically, given a choice of the intrinsic param- 
eters. For the BCV2 DTF, the intrinsic parameters are 
"00 7 "03/2 7 13, and /cut, and the extrinsic parameters are 
Ci 6 and the time of arrival tQ. 

In our tests, the maximization of the match is per- 
formed by way of a simplex-hascd search |20j in the con- 
tinuous space of intrinsic parameters (as opposed to the 
lattice-based searches used in bank-efficiency Monte Car- 
los); simplex methods have shown good efficiency in find- 
ing extrema in spaces of moderate dimensionality. In the 
light of the discussion of Fig. [2| in Sec. Ill Bl we under- 
stand that several values of the parameter B may yield 
local maxima of the match, corresponding to multiples 



Qp, and 2Qp of the target-signal precession fre- 
quency. To improve the robustness of our search, we 
run it repeatedly, starting with different initial values of 
B (100, 200, 300, 400, 500, and 600, covering the hkely 
range of B), so the final FF is usually picked out from the 
best of a few local maxima. In addition, after each run 
we restart the search from the current best-fit B, creating 
one more chance to escape a local maximum. (By con- 
trast, in Ref. [3 | we always started simplex-based max- 
imization at B = 20, too small compared to the values 
that we now see to correspond to the physical precession 
frequency. Because of this choice, Fig. 16 of Ref. 
shows more points clustered around what are likely to be 
values of B corresponding to f^p/2.) 

We test the performance of the BCV2 DTF for BH-BH 
target systems with component masses 



(7771,7772) 



12Mq 
IOMq 
SMq 

6M0 



> X 




(16) 



as well as for BH-NS systems with (7711,7712) 
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(IOMq, IAMq). We always take BHs to have maximal 
spins and NSs to be nonspinning. Without loss of gen- 
erality, we fix the directional parameters describing GW 
propagation and detector orientation, while we generate 
randomly 200 configurations of the other directional pa- 
rameters and of the local parameters of the binary (see 
Table I and the discussion around it in Ref. [31 for the 
definition of directional and local parameters). 

The FFs obtained are generally high, as seen by the 
plots of their distributions in Fig. |S1 That figure shows 
also the average fitting factor FF and the effective average 
fitting factor FFeff, as defined by Eq. (119) of Ref. [ij. 
The FFs have a strong dependence on the mass ratio 
of the target binary. The more asymmetric the masses, 
the hardest it is to get high FFs, probably because the 
number of precessional cycles is larger. Due to the im- 
proved search procedure, we find slightly better FFs for 
(10 + 1.4]Mo BH-NS binaries than the values obtained 
in Ref. [lj| (see Fig. llU|) . Moreover, we see that in Fig.[2| 
the relative number of points in the cluster correspond- 
ing to Op = rip is increased with respect to Fig. 16 of 
Ref. El. 

An important question is what ranges of DTF param- 
eters should be included in a template bank to be used in 
a search for a given class of target systems. Since all the 
BCV2 parameters are phenomenological, the straightfor- 
ward approach is to include templates corresponding to 
the range of the FF projection maps {i.e., to include the 
regions in the DTF parameter space where the maximum 
matches were achieved, for a representative set of tar- 
get systems). Figure El shows the (V'0,'03/2) and {ipo,B) 
sections of the FF projection. The clusters of maxima 
are labeled by the mass parameters of the target sys- 
tem; the spread in each cluster corresponds to different 
choices of the directional and local parameters of the bi- 
nary. The dotted rectangles show the regions used in Sec. 
IIV BI to estimate the number of BCV2 templates needed 
to search for the systems analyzed in this paper (with 
component masses (mi,m2) S [6, 12]M0 x [1,3]M0) and 
for heavier systems (with component masses (mi, 771,2) G 
[5, 2O]M0 X [5, 20]Mo). For comparison, the dashed poly- 
gon shows the {ipo, V'3/2) region suggested in Ref. [lj| for 
the heavier systems. 

Figure 13 shows the (V'Oj/cut) section of the projection 
map. The range of /cut has a weak dependence on the 
masses of the target system, very probably because the 
target signals end their evolution at the margin of the 
frequency band of good interferometer sensitivity (for 
LIGO-I), where little signal power is lost by excluding 
higher frequencies. Throughout this paper, we define the 
target signal ending frequency as the instantaneous GW 
frequency at the minimum energy circular orbit (MECO), 
as given by Eqs. (11)-(12) of Ref. [ll|, and plotted in 
Figs. 5 and 6 of Ref. |l|, and in Fig. 1 of Ref. [3. This 
frequency is inversely proportional to the total mass of 
the binary, and it is smaller for antialigned spin and or- 
bital angular momentum (i.e., negative n). 

In fact, it seems that the /cut parameter can be 



dropped altogether, as suggested by the following test: 
using the BCV2 DTF but fixing /cut = 400 Hz, we eval- 
uate FFs for (6-f 1)Mq and {12+3)Mq systems. The first 
mass configuration was chosen because, although it cor- 
responds to the largest ending frequency {i.e., the small- 
est total mass) among the systems studied in this paper, 
it shows the largest distribution of 7/13/2 in Fig. (SJ thus, 
by removing /cut from the maximization of the overlap, 
we could expect a very mild dependence on /cut, but 
a noticeable change in the FF projection. The second 
mass configuration was chosen because it corresponds 
to the smallest ending frequency {i.e, the largest total 
mass) among our systems. For both classes of systems 
the FF drops by only 0.5% with respect to searches that 
include /cut, with insignificant changes in the FF projec- 
tion ranges. 

Although the analytic maximization over the entire 
ranges of the extrinsic parameters carries little compu- 
tational burden, it is useful to constrain their ranges as 
tightly as possible to reduce the rate of false alarms. 
By constraining the ranges, we are in effect reducing 
the range of candidate signals that are compared against 
the detector output, and that have a (small) chance of 
being triggered by detector noise alone. In Fig. |S1 we 
show the (Ci,C2), {C3,Ci), and (Cs.Ce) section of the FF 
projection. The absolute magnitude of the Ck ranges is 
application-dependent, since it is determined by the over- 
all normalization of the waveforms.^ However, we notice 
that the C3,...^6 coefficients (which govern the amplitude 
of modulations) have magnitude comparable to the Ci^2 
coefficients (which multiply the unmodulated waveform) , 
and that the area occupied by the points shrinks slightly 
with decreasing total mass and more asymmetric mass 
ratios (corresponding to higher ending frequencies, and 
therefore greater signal power to be normalized). 

Asymmetric-mass-ratio binaries are not the only sys- 
tems having a large number of precession cycles. In 
Fig. 1^1 we show the number of precession cycles between 
40 Hz and the GW frequency of the test-mass ISCO, as 
evaluated from Eq. I^. In the upper panel, we consider 
asymmetric-mass-ratio binaries with m2 = IMq and mi 
ranging between IMQ and 15Mq; in the lower panel, we 
consider equal-mass binaries. We take two values of the 
opening angle; only the first object ("777,1") is spinning 
(maximally). We notice that the number of precession 
cycles is larger for binaries with asymmetric mass ratios 
than for equal-mass binaries, but that the largest number 
of precession cycles occurs for binaries with small, com- 
parable masses. For these, the BCV2 DTF has good per- 
formance: for instance, we find an average FF ~ 0.965 for 
(24-2)M0 systems (in the search, we fix /cut = 1000 Hz). 



^ Wc normalize all target signals and templates by assuming arbi- 
trary luminosity distances: this has no effect on the computation 
of FFs and template match metrics, which always involve nor- 
malized waveforms. 
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FIG. 6: FF projection maps onto the BCV2 intrinsic parameter space for target systems with component masses (mi, 7712) = 
{6, 8, 10, 12}AfQ X {1,2, 3}A'f0 and for (10 + 1.4)M0 BH-NS systems, with uniformly sampled directional and local angular 
parameters. The left panel shows the {'4>o,'4'3/2) section of the map, while the right panel shows the (11)0,3) section. The dots 
mark the values of the BCV2 parameters that achieve the maximum match for each target system. For each pair of masses, 
we draw an ellipse centered at the baricenter of the corresponding dot cloud, sized to include 90% of the dots; the axes of the 
ellipse follow the quadratic moments of the dots. The dotted rectangles show the regions used in Sec. IIV El to estimate the 
number of BCV2 templates needed to search for the systems with component masses (7711,7712) £ [6, 12]A/0 x [1,3]Mq and 
(mi, 7712) £ [5, 2O]M0 x [5, 2O]M0. In the left panel, the dashed lines enclose the region that was prescribed in Ref. [ij] for the 
heavier systems. 
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FIG. 7: FF projection maps onto the BCV2 (V'c/cut) pa- 
rameter subspace. See Fig. |H] for details. The cutoff seen at 
1000 Hz is due to a hard constraint arbitrarily imposed on 
the /cut search range. 



B. Performance of the BCV2P detection template 
family 



As we just saw, the BCV2 DTF offers good perfor- 
mance in matching the target waveforms. Its simple form 
allows the formulation of a simple prescription for tem- 
plate bank placement (as we shall see in Sec. IIV and 
the computational requirements are arguably economi- 
cal. However, it is not straightforward to extract physical 
information from the BCV2 DTF parameters, which are 
phenomenological. In this section, we discuss a modifi- 



cation of the BCV2 DTF for single-spin target systems, 
which is also written in the frequency domain, where the 
phenomenological intrinsic parameters t/^o, "03/21 £^nd B 
are replaced by the physical parameters mi, m2, X: and 

K. 

The first natural step toward a more physical 
parametrization is to replace the unmodulated phasing, 
■0o/"^/^ + V'3/2/^^/^ with the standard 2PN SPA phas- 
ing (as given, for instance, by Eq. (94) of Ref. ^^), which 
is a function of mi, m2, Xj ^^id k. The second step is 
to replace the power-law precession angle ap = Bf^^^^ 
with the analytic expression derived at Newtonian level 
Eq. ©, which is valid for single-spin systems, and which 
matches the numerical evolution of ap quite well [see 
Fig. (Q)] . In addition, since we expect this template fam- 
ily to be used in searches for asymmetric-mass-ratio bi- 
naries, for which the ending frequency falls at the margin 
of the band of good interferometer sensitivity, we do not 
include /cut among the intrinsic parameter, but instead 
we fix it to the GW frequency of the ISCO, as evaluated in 
the test mass limit. We denote this modified frequency- 
domain DTF as BCV2P (where the P in BCV2P stands 
for both "Pan" and "physical"). 

We emphasize that, although physically parametrized, 
the BCV2P DTF is not on the average much closer to the 
target waveforms than the BCV2 DTF. B oth D TFs suf- 
fer from the limitation emphasized in Sec. Ill CI namely, 
that the frequency components that appear in the mod- 
ulated GW phasing because of the evolution of the po- 
larization tensors do not simply occur at the precession 
frequency Hp, but also at its multiples (see Fig.^. Thus, 
the structure of the precession-frequency harmonics and 



11 



4 

2 



-2 

-4 
4 

2 ^ 
^ 

-2 

-4 
4 

2 



-2 

-4 



Ci3,5/10' 

-2 2 4-4 -2 2 4 -4 
12M0 T IOMq 









-2 2 4 -4 

SMq 






-2 2 4 

6M0 " 






4 

2 



-2 

-4 
4 

2 



-4 
4 

2 



-2 

-4 



-4 -2 2 4 -4 -2 2 4 -4 -2 2 4 -4 -2 2 4 

C,,3.5/10' 

FIG. 8: FF projection maps onto the BCV2 extrinsic parameter space, for the same target systems considered in Figs. |S|and 
Dots of different darkness mark the best-fit values of (Ci,C4), (€2,05), (CsjCe). The ellipses are drawn as in Fig. El 



the amplitude modulations arc not reproduced perfectly 
in either the BCV2 or the BCV2P DTFs. Truly exact 
physical templates for adiabatic spinning waveforms are 
so far available only in the time domain (where they are 
computed by solving the equations of motion) for single- 
spin binaries [T^ . 

We test the signal-matching performance of the 
BCV2P DTP only for {10+1 A)Mq BH-NS binaries. The 
BCV2P DTP is as effective as the BCV2 DTP, and much 
better than the standard unmodulated stationary-phase- 
approximated templates (SPA) . We plot the distributions 
of FPs for these three families in Pig.^| which shows also 
the average and effective average PP. In Pig.^] we plot 
FF as a function of the target-system parameter k. 

In Tab. we examine the parameter-estimation capa- 
bilities of M, rj, M., X K, by giving three character- 
istic quantities for each parameter: 

• the bias, defined as the average systematic error of 
the FF projection (it might be possible to remove 
the bias partially by careful characterization of the 
projection map); 

• the rms systematic error, caused by the spread in 
the PP projection due to the presence of unmodeled 
target-system parameters; 



parameter 


bias 


rms 


% within l-cr/3-O" 


M 


10.4% 


13.9% 


82%/98% 


V 


-0.104 (abs) 


0.187 (abs) 


78%/98% 


M 


1.2% 


1.3% 


60%/100% 


X 


-20.0% 


15.1% 


80%/96% 




0.020 (abs) 


0.153 (abs) 


80%/98% 



TABLE I: Bias, systematic rms error, and percentage of es- 
timators falling in the 1-a and 3-a intervals for the BCV2P 
DTF parameters M, rj, Ai = Mrf^'^ (the chirp mass), x and 
K. The rms errors for M, M, and x sue given as percentages 
of the target-system value of those parameters. 



• the Gaussianity of the distribution, as character- 
ized by the percentage of estimators falling in the 
1-a and S-cr intervals (for a Gaussian distribution, 
these should be ~ 69 % and ^ 100 %, respectively). 

These systematic errors are distinct from statistical error 
due to detection in noise, which is roughly inversely pro- 
portional to source strength (i.e., S/N). We see that chirp 
mass is the parameter that can be estimated most pre- 
cisely, with ^ 1% bias and rms deviation. Interestingly, 
K can also be estimated rather well. 

Despite its fine parameter-estimation performance, the 
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FIG. 9: Number of precession cycles for asymmetric-mass- 
ratio binaries (upper panel) with m2 — IMq, and for equal- 
mass binaries (lower panel), as functions of mi, for two values 
of the opening angle k. 
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FIG. 10: Distribution of BCV2, BCV2P, and 2PN SPA fitting 
factors for (10-f 1.4)M0 BH-NS target systems with uniformly 
sampled directional and local angular parameters (400 sets). 
See Fig. |S] for details, but note that here FF bins are shown 
on a linear scale. 



BCV2P DTF has some disadvantages with respect to the 
BCV2 DTF and to the time-domain physical templates 
of Ref. For example, it will be more complex to build 
the template match metric and to place down templates 
for BCV2P than for BCV2 (see Sec. ITvH below). From 
the point of view of FF, false-alarm rate, and parameter 
estimation, the BCV2P DTF is less attractive than the 
physical templates of Ref. However, these might well 
be too computationally burdensome to be implemented 
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FIG. 11: Average fitting factors achieved by the BCV2, 
BCV2P, and 2PN SPA DTFs for (10 + 1.4)M0 BH- 
NS target systems with uniformly sampled directional and 
local angular parameters (400 sets), plotted against k. 
The vertices of the segmented curves show the FFs av- 
eraged on the sets of target systems with k in the bins 
[-1.0, -0.8), [-0.8, -0.6), . . . , [0.8, 1.0]. The bars show the 
sampling error on the bin averages. 



for online searches; in that case, the BCV2P DTF could 
be use as an efficient first stage in a hierarchical search 
strategy. 



IV. A PROCEDURE FOR TEMPLATE 
PLACEMENT USING THE TEMPLATE MATCH 
METRIC 

In this section we show how to place the BCV2 tem- 
plates within a certain DTF parameter re gion, w hile 
guaranteeing a chosen minimum match MM [2^ |23 de- 
fined by 



MM = min max p(A, A') 

A A'Gbank 



mm max 



(MA), MA')) 



:(17) 



A A'Gbank ^ {h{X) , h{\)) {h{\') , h{X')) 

(where p is the match) in terms of the noise inner product 



Jo 



Snif) 



(18) 



with Snif) one-sided noise spectral density (given 
for this paper by Eq. (68) of Ref. 0). Although the 
maximization over the extrinsic DTF parameters can be 
carried out analytically, the existence of extrinsic param- 
eters still influences the placement of templates, as dis- 
cussed, e.g., in Sec. VI of Ref. 0|. 



A. Template metric of the BCV2 DTF 

The match between templates with close parameter 
values can be approximated using a metric in parameter 
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space [11111113, 

p{\ A + AA) = 1 - .gcijAA^AA^ 



(19) 



The components of the metric can be expressed in terms 
of first derivatives of the template waveforms. 



9CD 



1 I dh dh 

1 I dh 



{h\h)^\dX<^' 



^'aA^/ 



(20) 



It should be noted that the match between nearby tem- 
plates with different cutoff frequencies cannot be de- 
scribed by the metric of Eq. H2U|) : in this paper, we shall 
consider only the problem of placing templates that share 
the same cutoff frequency /cut- For binaries with low 
masses, the waveforms end at relatively high frequencies 
compared to the band of good interferometer sensitivity, 
so /cut does not play an important role. (For instance, 
adopting the language of Sec. VI of Ref. 0, /cut can 
be set by requesting that (/i(/cut), h{+oo)) ~ 0.99, which 
yields /cut — 400 Hz using the Newtonian amplitude evo- 
lution /-■^/^) 

To evaluate the metric gcD, we project 
/i(V'o, ■03/2, -S, /cut, Cfe;/) [Eq. ijH)] onto an orthonor- 
mal basis ei{f) {j = 1, . . . , 6), writing 

h = A'e,{f), (e„e,) = %; (21) 

the template is then normalized if 

{h,h) ^ AjA^ ^1. (22) 

A convenient choice of the basis functions £.;(/) is the 
following: obtain the functions e'i.3,5(/) from the Schmidt 
orthonormalization procedure. 









an 




1 








63 




031 033 




cos(6/ 


-2/3) 


/-^/« (23) 




. ^5 . 




151 053 055 




sin(S/- 


-2/3) 


1 



(note that the are functions of B); define 62,4,6 (/) 
from 

e„+i(/) = i ?„(/), 71 = 1,3,5; (24) 
(for / > 0); finally, define 

e„(/)=^n(/)6'*(/)+2"^*=. (25) 

Our parameter set is now A*^ = {A*;.t"} = 
{A*; tc, B, tpo, 1P3/2}, with i = 1, . . . , 6 and a = 0, 1, 2, 3. 
The A\ along with x° = tc, are extrinsic parameters, 
while the (for a = 1, 2, 3), are intrinsic parameters. 
Equation (|22|l is complemented by the useful relations 



/ dh dh 

XdA^'dAl. 
I dh dh 

X'dA^'dx^. 
I dh dh 

\dx°^''dxP/ 



A^, 

0, 



aaI— — 



\ = 



(26) 
(27) 
(28) 
(29) 
(30) 



also, the first derivatives of 6„(/) with respect to the 
can be summarized in the differential expression 



de„ = 27ri/e„ dtc + if ^^^e„ dipo 
+*/-2/3e„d03/2 + ^|^dS 



„i*(/) + 277i/t. 



(31) 



Assuming AiA^ = 1, we write the match between nearby 
templates as 



piX,X + AX) = 1-5c'd(^\S)AA^AA^ = 1 



1 r 
2 



AA' Ax" 



In language of Ref. ^1|, gcD is the full metric, which 
describes the match between nearby templates in terms 
of differences between all their parameters. In general, 
gcD can depend on all intrinsic and extrinsic parameters, 
but in our case it depends only on A' and B, but not on 
■00,3/2 s-ncl tc [as a consequence of Eqs. 



Sij AiAj 






^'^'"\dx'^' dxf'l - 



AA^ 
Ax^ 



(32) 



r 



To determine the spacing of bank templates along the 
intrinsic-parameter directions, we work in terms of the 
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extrinsic-paramcter-maximizcd match 
= max p^,^{A\x^;A^+AA\t, + At„x"' + Ax''). 



(33) 



which is approximated by the projected metric 5?™^ , cal- 
culated by maximizing Eq. H32|l over AA^ and Ate, while 
fixing Aa;" (see Eq. (65) of Ref. 0|). Maximizing first 
over AA^ , we notice that the submatrix Sij — AiAj is 
degenerate, with the single null eigenvector A^ [this de- 
generacy occurs because the match (|32|) remains constant 
when AA^ moves parallel to A^ — it must then be broken 
when we impose, without loss of generality, AkAA'' = 0]. 
We find the maximum 



l^^[A'A^G,ja.p{B)] Ax" Ax" 



at the location 



AA" 



(34) 



(35) 



[note that AkAA'" = due to Eq. The tensor G^jaf3 

of Eq. (|34|l is given by 



de 



Further maximizing H34(l over At^, we find the maximum 



Pmax(A^ x^; a;^ + Ax^) = 1 - glfiA^,B) Ax^'Ax'^ , 



(37) 



where g^f^^ is the three-dimensional projected metric 



A'A^G 



A'A^A^A'^Co&G 



ij&P 



A'A^G 



ijOO 



In general, the projected metric can depend on all the 
parameters, but in our case it depends only on A* and B. 

We can now go back to the problem of template place- 
ment. Within the metric approximation, if we choose a 
target template [the h{X) of Eq. p7|l ] with parameters 
A'-^ = {A^,x'^), the nearby templates with (extrinsic- 
parameter-maximized) match greater than MM must 
have intrinsic parameters x'^ + Ax^ that lie within the 
ellipse £{A^,B, MM) specified by 

£{A3 , B, MM) : gPjj {A^ , B) Ax^Ax^ < 1 - MM . (38) 

The shape of this ellipse depends on the target-template 
extrinsic parameters A^ , which is not appropriate for the 
placement procedure that we are seeking, which should 
be formulated in terms of the intrinsic parameters only. 

To that purpose, Ref. suggested adopting the aver- 
age shape of the match ellipse, as obtained by averaging 




FIG. 12: Strategies to place templates in the space of BCV2 
intrinsic parameters. 

Top. Suppose we have already placed a template at xf, and 
consider only two sets of target-template extrinsic parame- 
ters, associated with which are two different match ellipses, 
£i and £2- In order for the bank to guarantee the MM for 
both sets of extrinsic parameters, each of £1 and £2 must con- 
tain at least another template at a location different from . 
This can be achieved with two templates (X3 and X4 in the 
figure), or more easily (but less optimally) with a single tem- 
plate {x2 in the figure) in the intersection of £1 and £2- 
Bottom. In the idealized situation where all ellipses have the 
same elongated shape, but take all possible orientations, plac- 
ing templates along a multilattice (black dots) can be much 
more efficient than tiling on the basis of the minmax region 
(gray). To obtain the multilattice, we construct a set of 
maximized-match ellipses evenly separated by angles 8 ~ b/a, 
so that the circle (dashed) through the intersections of the el- 
lipses has radius ~ a. The unit cell of the multilattice is then 
given by the intersections of the semimajor axes of the ellipses 
with the square inscribed in the dashed circle. 



the maximized match [and hence the left-hand side of 
Eq. H38() ] over the extrinsic parameters. In Ref. the 
averaging weights were determined from the prior dis- 
tribution of the target-signal physical extrinsic parame- 
ters, in such a way that a template placement procedure 
guided by the average match contour would guarantee a 
certain expected detection efficiency. 

In our case, however, the extrinsic parameters are not 
physical, and do not have obvious prior distributions. 
We take a conservative approach, and we require that 
for every value of the intrinsic parameters x'^ in the bank 
and for every possible value of the extrinsic parameters 
A^ , there exist a nearby bank template with Ax'^ within 
the ellipse £{A^ , B,MM). In principle, the neighboring 
template that satisfies the criterion could be different for 
different A^ , as illustrated at the top of Fig. ^] For 
simplicity, however, we adopt a suboptimal strategy, re- 
quiring the existence of a nearby template within the 
intersection of all the ellipses centered at x'^ (the gray 
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region at the top of Fig. I12() . 



^mmmax 



(S,MM) = f]£iA\B,MM). 



(39) 



A3 



In analogy to the usage of Ref. |2j|, we denote 
fminmax('S, MM) as the minmax region, because it cor- 
responds to considering the contours of the match max- 
imized over the extrinsic parameters of search tem- 
plates, and minimized over the extrinsic parameters of 
the prospective target signals, 



mm 

A3 



> MM. 



(40) 



The minmax regions £'minmax(S, MM) may no longer be 
elliptical in shape, but their linear dimensions still scale 
linearly with Vl — MM, as it happened for the ellipses 
£{A^,B, MM). 

It is worth pointing out that if the orientations of 
the ellipses £{A^ ,B,MM) vary substantially when A^ is 
changed, relying on different nearby templates to achieve 
the MM for different extrinsic parameters can be much 
more efficient than relying on the minmax region. As 
a simple example, consider a situation in which all the 
ellipses £{A^ , B,MM) have the same shape, with semi- 
major and semi-minor axes a and b (with a ^ b), but 
assume all possible orientations, as shown at the bottom 
of Fig. El In this case, the minmax region has area ^ b^. 
If we place templates according to the minmax prescrip- 
tion, the template density becomes ~ 1/6^, much larger 
than the density 1/ (ab) associated with each individual 
ellipse. 

On the other hand, we could place templates on the 
sides of squares with size ~ a, separated by a param- 
eter distance ^ b, as shown at the bottom of Fig. [T^ 
In this multilattice, the average area occupied by each 
template is ^ a^/{a/b) ^ ab, corresponding to a density 
~ 1/ (ab) , much better than obtained with the minmax 
prescription. Unfortunately, generating the appropriate 
multilattice is definitely more complicated than using the 
minmax prescription; it is also not clear whether for the 
BCV2 DTF we have in fact elongated ellipses with dra- 
matically different orientations. Thus, we will adopt the 
minmax prescription in the rest of this paper. 

We can approximate the minmax regions for the BCV2 
DTF by sampling the A^ randomly. In Fig. ^1 we plot 
an example minmax region for B = 60, /cut = 400 Hz. 
This figure is typical in the sense that the minmax re- 
gions in the {x^,x'^,x^) = (V'o, '03/2) 'S) space can be ap- 
proximated rather well by ellipses. As a consequence, 
we can rely on (yet another) metric g^fj{B) in the space 
of intrinsic parameters, whose match ellipses lie within 
the corresponding minmax regions, but have similar vol- 
umes. Once we are equipped with gg^p{B), the placement 
of the BCV2 templates can be performed along the lines 
of conventional template placement procedures. 




-100 Ay. 



3/2 



FIG. 13: The BCV2 DTF minmax region for B = 60, /cut = 
400 Hz (and any ipo, "03/2), obtained by sampling Eq. H40|l 
with a random distribution of . 



AB 




in the BCV2 intrinsic- 
The orthonormal basis 



FIG. 14: Template placement 
parameter space {■tpo,i>3/2,B). 
('/(I) 1 '7(2) ' V{3) ) is constructed in such a way that 77"^) and 77^2) 
both lie in the (00, '03/2) plane, along which the projected 
metric is constant. Unit cells constructed along these basis 
vectors, as shown in the figure, can be extended to the en- 
tire intrinsic-parameter space consistently. The coefficient e 
is v'4(l -MM)/3. 



B. Template placement 

The standard local prescription for template placement 
is to follow a cubic lattice"' constructed with three or- 
thonormal basis vectors ^ V"2y V^^)} 1 setting the side 
length [as measured with the metric gg^^iB)] of the unit 



* Under some circumstances, other lattices can provide better 
packing: for the BCV2 DTF, preliminary tests suggest that a 
tctrahcdral lattice could reduce the number of templates by one 
fourth with respect to a cubic lattice. 
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FIG. 15: Effective parameter volume of a single template cell 
as a function of B, assuming a cubic lattice with MM = 0.97. 
We fix /cut = 400 Hz. 



systems with component masses (mi,TO2) € [6, 12]AfQ x 
[1, 3]Mq, we select (following Fig.© the BCV2 DTF pa- 
rameter region 



105 



•03/2 
103 



B 

102 



G [1.5, 8.5] X [-4.5, 0.5] x [0, 10] ; (43) 



our estimate for the nmirber of templates is A/'bh-ns — 7 x 
10^ for MM=0.97. By contrast, for BH-BH systems with 
component masses {mi, 7712) £ [5, 20]Mo x [5,2O]Af0, we 
select the BCV2 DTF parameter region 



tpo V'3/2 
IQS' 103 



B 

102 



e [0,2.5] X [-2.5,1] X [0,4.5] ; (44) 



out estimate for the number of templates is A/'bh-bh — 
8 X 10"* for MM=0.97. However, it is not clear that 
precessing-binary templates are needed for this entire re- 
gion: for binaries with relatively high component masses, 
"BCVl" (unmodulated) templates can already yield high 
FFs. 



cell equal to i/4(l - MM)/3 [Tl]. For general metrics 
that depend on location in parameter space, such a local 
lattice cannot usually be extended consistently to cover 
the entire space. Luckily, this is possible in our case be- 
cause of the translational invariance of g^p{B) along the 
■00 and ■03/2 directions. 

We first identify a set of orthogonal basis vectors 
{^(1)' '^p)' ^(3)} each point x", with the property that 
both 77"^^ and r]"^^ lie within the (^O)'03/2) plane. One 
such set follows from defining 




(41) 



this Ty^gj is orthogonal to all tangent vectors that lie 
within the ('0o,'03/2) plane. We can complete the ba- 
sis with any pair of ^-orthonormal vectors I??"!): '7(2)} 
that plane. Due to translational invariance, the cubic 
lattice constructed with this basis can be extended con- 
sistently along the ipo and ^03/2 directions, covering a thin 
slice of parameter space parallel to the ("00, "03/2) plane, 
with coordinate thickness 



AB 



4(1 -MM) 3 

5 '^(3) 



4(1 - MM) 



(42) 



along the B direction. We can stack these slices to cover 
the entire tridimensional parameter space (see Fig. I14|l . 
In Fig. 1151 we plot the parameter volume of the cube 
inscribed in the MM=0.97 minmax regions {i.e., the ef- 
fective volume around each template in the cubic lattice) 
as a function of B. 

We are now ready to give a rough estimate of the num- 
ber of templates required for a matched-filtering search 
for BH-NS and BH-BH target systems. We assume a 
single value (400 Hz) for /cut- For BH-BH and BH-NS 



V. CONCLUSIONS 

Reference Q introduced the BCV2 DTF for use in 
precessing-binary searches; the BCV2 templates are writ- 
ten directly in the frequency domain, and depend on 
ten phenomenological parameters (the intrinsic param- 
eters ipQ, 03/2, B, and /cut, and the extrinsic parame- 
ters Ci^,,,_e). The modulational effects of precession in 
the target systems are modeled by the Ci^..._6 and by 
the single intrinsic parameter B; this marks a definite 
improvement over previous attempts at spinning-binary 
templates 0, S H Ellll El [13, which require sev- 
eral intrinsic parameters to model precession, leading to 
larger and more computationally expensive search tem- 
plate banks. 

In this paper, we have tested the signal-matching per- 
formance of the BCV2 DTF against target signals from 
asymmetric-mass-ratio binaries with component masses 
(7711,7772) G [6, 12]M0 X [1,3]M0 (with both components 
maximally spinning), and from (10-1- 1.4)AfQ BH-NS bi- 
naries (with maximally spinning BHs and nonspinning 
NSs). The waveforms were computed at the 2PN order 
in the adiabatic approximation [j] . We found very good 
fitting factors, averaging between 0.94 and 0.98; our re- 
sults are summarized in Fig. \S\ By means of the FF pro- 
jection map (see Sec. IIII A|) we also identified the region 
in BCV2 DTF parameter space that must be included 
in a template-bank-based search for these systems; our 
results arc shown in Figs. Inland 

Because all the asymmetric-mass-ratio binaries consid- 
ered in this paper generate waveforms with ending fre- 
quencies {i.e., MECO frequencies, see Eqs. (11)-(12) of 
Ref. [3) at the margin of the band of good interferom- 
eter sensitivity, we are free to fix /cut to a reasonable 
value (see, for instance, the discussion below Eq. ^U^), 
reducing the dimensionality of template space without a 
corresponding degradation in the FFs. 
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In addition, by a closer study of precessional dynamics 
and GW generation in single-spin binaries, we were able 
to relate the BCV2 phenomenological parameter B to the 
physical parameters of the target binary (see Sec. Ill B)l . 
In the process of doing so, we realized that waveform 
modulations occur at the fundamental and at higher and 
lower harmonics of the precession frequency fJp (for our 
target systems, the frequency at which the orbital an- 
gular momentum and the spin precess around the total 
angular momentum; for the BCV2 templates, the equiva- 
lent precession frequency corresponding to a choice of B). 
In the target signals, the higher harmonics arise because 
modulation is caused by oscillations in the components 
of the polarization tensor [see Eq. l(T^ ]. and not directly 
by the precession of the orbital angular momentum and 
spins. In the BCV2 templates, the higher harmonics fol- 
low naturally from Eq. Q, and can even be made dom- 
inant by an appropriate choice of the phenomenological 
coefScients Ci^..._6 (see Fig. E)). 

These considerations allowed us to understand cer- 
tain features in the distribution of best-fit B against the 
target-system opening angle k that had remained unex- 
plained in Ref . Il4l | . The analysis performed in this paper 
suggests also a modification of the BCV2 DTF, whereby 
the three phenomenological parameters "00 j "03/21 f-nd 
B are replaced by the four physical parameters M, 77, 
K and X- This modified DTF (BCV2P) has a signal- 
matching performance comparable to (or slightly better 
than) the BCV2 DTF, for (10 + IA)Mq BH-NS bina- 
ries (see Figs. llUland lllll : we expect that we would find 
slightly higher FFs if we were to include higher-order PN 
corrections in ap. The BCV2P DTF has the advantage of 
providing the straightforward and reliable estimation of 
certain parameters of the target system, such as M and 
K (see Table nj; it has the drawback of depending on four 
rather than three intrinsic parameters, and its template 
match metric would be even more complicated than the 
BCV2 metric described in Sec. IIV Al Therefore, we sug- 
gest that the BCV2P DTF could be used for follow-up 
searches on a reduced set of BCV2 triggers, or to estimate 
target-system parameters, if it turns out that we cannot 
computationally afford the physical template family of 
Ref. which is also written in terms of the physical 
parameters of the target system, and which yields FFs 



very close to unity. 

Last, we computed the full ten-dimensional template 
match metric in the (0Oj "03/21 ^c, B, Ci, . . . ,Cg) param- 
eter space, and the tridimensional projected metric in 
the (00, 03/2, B) subspace, obtained by projecting out 
the seven extrinsic parameters. The projected metric 
does not depend on tpo and 03/2, but only on B (and 
/cut)- We described a prescription to place BCV2 tem- 
plates in a search bank using minmax regions^ and ex- 
ploiting the (00, "03/2) translation invariance of the met- 
ric. Fixing /cut = 400 Hz and adopting a cubic lattice, 
we find that the estimated number of templates required 
for MM = 0.97 is ~ 7 x 10^ for target systems with 
component masses (7711,7712) € [6, 12]Mo x [1,3]Mq, and 
8 X 10"* for target systems with component masses 
(7771,777.2) e [5, 20]Afo X [5, 2O]M0. These numbers make 
it unlikely that a straightforward BCV2 search could be 
performed online^ with the current computational re- 
sources, for these full parameter ranges. Workarounds 
might include carefully isolating (by extensive Monte 
Carlo runs) the subregions of parameter space where 
spin effects are weaker, so that nonspinning templates 
(such as the BCVl DTF |||) can be substituted effec- 
tively for BCV2 templates, or employing the BCVl and 
BCV2 DTFs sequentially in a hierarchical search. 
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